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ABSTRACT 



We present an analysis of super-giant molecular complexes (SGMCs) in the overlap region of the Antennae galaxy merger, based on 
ALMA CO(3-2) interferometry and VLT/SINFONI imaging spectroscopy of IT 1-0 S(l) at angular resolutions of 0.9" and 0.7", 
respectively. All but one SGMC have multiple velocity components offset from each other by up to 150 km s _1 . H 2 line emission 
is found in all SGMCs and the kinematics of IT and CO are well matched. IT/CO line ratios vary by up to a factor of 10 among 
SGMCs and different velocity components of the same SGMCs. We also identify the CO counterpart of a bright, compact source 
of near-IR IT line emission, which shows no Bry, and was first identified with SINFONI. This source has the highest IT/CO line 
ratio, and coincides with the steepest CO velocity gradient of the entire overlap region. With a size of 50 pc and a virial mass of a 
few 10 7 M B it is perhaps a pre-cluster cloud that has not yet formed significant numbers of massive stars. We present observational 
evidence that the H 2 emission is powered by shocks, and demonstrate how the IT 1-0 S(l) and the CO(3-2) lines can be used as 
tracers of energy dissipation and gas mass, respectively. The variations in the IT/CO line ratio may indicate that the SGMCs are 
dissipating their turbulent kinetic energy at different rates. The compact source could represent a short (~ 1 Myr) evolutionary stage 
in the early formation of super-star clusters. 
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1. Introduction 

Major gas-rich mergers are important sites of star formation and 
galaxy evolution in the Universe. The Antennae galaxy merger 
(NGC4038/4039) is an ideal target for studying in detail how 
galaxy interactions affect the interstellar medium and star for- 
mation. Most stars in the Antennae form in sup er-star clusters 
(SSC s) with stellar masses up to a few 10 6 M Q tiWhitmor e et alj 
l2010l) located where the two galaxies permeate each other, the 
'overlap region'. Super-giant molecular complexes (SGMCs) 
with masses of several 10 8 M G and sizes of ~500 pc have been 
identified i n CO( 1 -0) in the overlap region with th e OVRO inter- 
ferometer (I Wilson et al.ll2000l) . IUeda et al.l d2012l) have recently 
reported higher resolution (~ 100 pc), CO(3-2) observations of 
the Antennae obtained with the SMA. 

The formation of SSCs involves a complex interplay of 
merger-driven gas dynamics, turbulence fed by the galaxy in- 
teraction, and dissipation of the kinetic energy of the gas. 
Hydrodynamic simulations suggest that massive complexes of 
cold gas, akin to SGMCs, form where gas flow s trigger com- 
press ion, cooling and gravitational fragmentation (Tevssie r et alj 
2010). Within SGMCs, a hierarchy of structures must form 
including clouds tha t are sufficiently massive to form SSCs 
(IWeidner et alJl20lOh . 

Recent VLT/SINFONI imaging spectroscopy of the peak of 
pure-rotational H2 e mission in the overl ap region previously ob- 
served with Spitzer (He rrera et alj|201 H HI 1 hereafter) revealed 
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bright diffuse H2 line emission associated with an SGMC and a 
compact (~0.6", ~50 pc) source. HI 1 proposed that the compact 
source may be a massive cloud on its way to form a SSC within 
the next few Myr. The H2 lines are powered by shocks and trace 
energy that is being dissipated and radiated away as the cloud 
complex, and a pre-cluster cloud (PCC) within, grow through 

gas accre tion. 

iHerrera et al. l d20Tl had sub-arcsecond resolution SINFONI 
data of shocked gas, but lacked CO observations at similar res- 
olution, which are required to probe the bulk of the gas on the 
relevant scales of <100 pc. In this letter, we take advantage of 
the recently released ALMA science verification observations of 
CO(3-2) in band 7 (-345 GHz) at ~1" angular resolution to 
compare the morphology and kinematics of CO(3-2) and H2 
line emission in the Antennae overlap region. We show that to- 
gether the ALMA and VLT observations provide the comple- 
mentary (mass and energetics) data needed to characterize the 
dynamical state of SGMCs and to search for pre-cluster clouds. 



2. Comparison of ALMA and SINFONI observations 

Our analysis relies on two data sets. First, CO(3-2) line emis- 
sion was obtained during ALMA science verification. These data 
are part of a mosaic of the Antennae obtained in 10 hr of ob- 
serving time in band 7 (345 GHz) between May and June 201 1, 
with 10 to 13 antennae and baselines from 25 to 200 m. This 
gives an angular resolution of 076xl''l (66 pcxl 15 pc, at a dis- 
tance of 22 Mpc) and covers the entire overlap region. The data 
have an intrinsic spectral resolution of 0.85 km s and were 
binned into channels of 10 km s _1 . We used the reduced data 
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R.A. (J2000) R.A. offset [orcsec] from 12:01:55.12 (J2000) 

Fig. 1. Left. ALMA CO morphology shown on top of our CFHT K-band continuum image (HI 1). Dotted boxes mark the two ALMA 
mosaics, the solid box marks the overlap region. Right. H2 1—0 S(l) morphology as seen with SINFONI. Boxes mark individual 
SINFONI fields-of-view, contours show CO(3-2) from 2 to 42 Jy km s~' beam -1 in steps of 8 Jy km s~' beam -1 . The inset at 
the bottom left of the right panel shows the ALMA beam. We also mark massive and young SSCs (asterisks), and the compact H2 
source PCC discussed in §[3] 



cubes publicly available on the ALMA websit^], which we cor- 
rected for the primary beam attenuation. Low spatial frequencies 
were filtered out because of missing short spacings that cause a 
loss of extended structures >4"and negative sidelobes adjacent 
to bright emission. To measure fluxes we used a clipped cube 
where all pixel values <2rx (6 mJy beam -1 ) were set to zero. To 
quantify the missing flux, we constructed spectra at the center of 
the overlap region with the angular resoluti on of the single-dish 
JCMT and HHT CO (3-2) observations of lZhu etalJ d2003l) and 
ISchulz et alJ (12007), respectively. We found about half the total 
flux of the single-dish data. The total flux of the SG MCs in the 
overl ap region agrees with the SMA observations dUeda et al. 
120121) . Second, we used VLT/SINFONI imaging spectroscopy 
of ro-vibrational H2 lines in the near-infrared K-bwd at R=3000 
in four regions each 8"x8" in size (Fig.Q]). They were observed 
in February 201 1 with on-source integration times of 40 minutes 
per pointing. Our previous observations of SGMC 2 have been 
discussed by HI 1 . We obtained and reduced the data of the three 
additional fields in a similar way. 

Fig. Q] shows the spatial distribution of the H2 and CO emis- 
sion in the Antennae overlap region at a spatial resolution of 
<1" (100 pc). The right panel displays the color image of the 
H2 1-0 S(l) emission with CO(3-2) contours for comparison. 
Intensity maps were obtained from the SINFONI and ALMA 
data cubes by fitting Gaussian profiles to the spectra in each spa- 
tial pixel, using one component for SINFONI and up to three 
components for the higher spectral resolution ALMA data. In 
Fig. [2] we compare the velocity fields of these two lines. For the 
CO velocities we computed the first moment map, and for H2 we 
constructed the velocity map from Gaussian fits. 

The ALM A channel maps a re very similar to the SMA ob- 
servations in lUeda et all (120121) . We identified all SGMCs dis- 
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covered by IWilson et all d2000l) except for SGMC 3, which is 
not covered by ALMA. The new data show two velocity com- 
ponents in SGMC 1 and 2. We also decomposed SGMC 4/5 into 
two smaller complexes at positions (-3,-9.5) and (-3,-3.5) in 
Fig- E eacn °f which has two velocity components. We kept the 
name SGMC4/5 for the first complex, and labeled the northern 
extension SGMC 6. CO and H 2 spectra of each SGMC are dis- 
played in Fig. [2] The CO spectra were integrated over each box 
in Fig. Q]using the clipped cube. 

Table[T]lists the CO line properties of each SGMC and each 
velocity component, named a and b. Fluxes, velocities and line 
widths of all components are measured from Gaussian fits to the 
spectra in Fig. [2] Error bars include only fit uncertainties, not 
the systematic errors owing to the missing short spacings. We 
estimated the R3-2/1-0 = ^co(3-2)/^co(i-0) ratios of the SGMCs 
by comparing the ALMA and OVRO data, after smoothing to a 
common resolution. Iqo is the integrated intensity in K km s . 
This ratio varies from source to source be tween 0.3—0.8, with 
a mean value of 0.5, consistent with what Schu lz et al. I d2007h 
found from single-dish data for the entire overlap region, as well 
as the peak line ratios measur ed with the SMA for each individ- 
ual SGMC (lUeda et alj|2012h . 

We estimated gas masses from the CO fluxes, where the Xco 
factor is the main source of uncertainty. To be consistent with 
previous studies, w e used the same Xco factor for CO(1-0) as 
IWilson et all d2000), X co =3 x 10 20 H 2 cm' 2 (Kkm s~'r ', and 
adopted the scaling R3_2/i-o=0.59 of ISchulz et alJ d2007l). Our 
mass estimates (Tab. [TJ are comparable to those of IWilson et al.l 
d2000h . Similar to Hll, our data show that most of the H 2 
emission away from SSCs is powered by shocks, not UV heat- 
ing in PDRs. The observed H 2 2-1/1-0 S(l) ratios, 0.1-0.2 
(Tab. [T]) can be accounted for by PDR models, but only for 
high UV fields (y>10 4 , in units of the mean val ue in the solar 
neighborhood) and high densities («h>10 5 cm -3 , Le Petit et al. 
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Fig. 2. Left panel: Velocity map of H2 1-0 S(l) (left) and first-moment map of CO(3-2) (right). Dotted lines mark the position- 
velocity cut shown in the right panel. Mid panel: Integrated line profiles of CO(3-2) (red) and H2 1-0 S(l) (black) for each SGMC. 
Blue spectra show the CO(3-2) lines convolved to the spectral resolution of SINFONI. Black bars correspond to 5 xl0~ 16 erg s _1 
cm -2 A -1 for H2 1-0 S(l) and red bars to 1.5 Jy for CO(3-2), for all SGMCs. Offsets between the black and blue spectra indicate 
variations in the CO-to-Hb ratio between velocity components. Right panel: H2 1-0 S(l) position-velocity diagram of SGMC 2 
(HI 1), with CO(3-2) emission shown as contours in steps of 0.02 Jy beam -1 starting at 0.03 Jy beam -1 . 



2006). These conditions exist near massive stars embedded in 
molecular clouds, but the extended H2 emission is generally 
not observed to peak near the brightest SSCs. The notable ex- 
ception is the embedded cluster in SG MC4/5, which, however, 
has a higher line ratio dGilbert et al.l 120001) . In addition, the 
[CII]157jL/m/[OI]63jum line ratio indicates that the mean radi - 
ation field in the overlap region is ^~10 3 (ISchulz et al.l 12.007b . 
an order of magnitude smaller than that required to account for 
the H2 2-1/1-0 S(l) ratios. The other luminous SSCs have al- 
ready dispersed most of their gas and dust. Thus, the brightest 
UV-heated gas would come from clouds away from the clusters. 
Outside the clusters the intensity of the radiation field is not very 
high. For a 5xl0 6 M SSC, the UV field is *>10 4 only out to 
a distance of <100 pc from the cluster and less if we include 
extinction. 

Our analysis has three main results. (1) All SGMCs have 
H2 1-0 S(l) line emission and the H2 1-0 S(l ) kinematics 
match those of CO(3-2) well. IZhu et all (120031) and lSchulz etail 
(120071) found that the CO is emitted from gas at temperatures 
with T~30-150 K. The excitation temperature of the warmer 
H 2 gas emitting in the near-IR is -1000-2000 K (Hll). The 
similarity of the gas kinematics of CO and H2 indicates that 
warm and cold gas are closely associated. (2) Figs.Q]and[2]show 
large variations in the H2/CO ratio between SGMCs and between 
individual velocity components in the same SGMC (Tab. Q]). 
Extinction is an unlikely cause, because the near- and mid-IR 
H2 pure-rotational emission-line regions have similar morpholo- 
gie s (at least at the 5" scales resolved by Spitzer-IRS; see Fig. 4 
in [Brandl et af] |2009l and the discussion in Hll). In Hll, we 
related the H2 emission to the dissipation of kinetic energy. With 
this interpretation, the H2/CO ratio traces differences in the en- 
ergy dissipation rate per unit mass, which must be related to the 
dynamical evolution of the gas. (3) All clouds have two spatially 
separated velocity components as seen in the channel maps of 
SGMC 2 in Fig. [3] In the other clouds, the spatial offsets be- 
tween velocity components are less obvious, possibly because 



of projection effects. The velocity difference between compo- 
nents within an individual SGMC is up to 150 km s -1 (Fig. [2J. 
Given the size and mass of the SGMCs, this is too large to be 
accounted for by the gas self-gravity. The gas kinematics is most 
likely driven by the galaxy interaction. Single-dish observations 
found similar velocity gradients in the extended emission around 
the SGMCs, which further supports this idea. The two compo- 
nents in SGMC 2 are spatially separated and resolved by ALMA. 
We estimated their sizes and found their virial masses to be com- 
parable to the molecular masses derived from the CO fluxes. 

3. The compact H 2 source 

The ALMA maps also give new insight into the nature of the 
bright, compact H2 emitter associated with SGMC 2, PCC, re- 
cently discovered by Hll. PCC is the brightest H2 line emitter 
in the overlap region. It is not detected in the lower-resolution 
CO(1-0) data, but is an emission peak in the SMA and ALMA 
maps (Fig. [TJ. Isolating the PCC CO counterpart from the sur- 
rounding extended emission is difficult with an algorithm like 
CLUMPFIND because of the large velocity gradient across PCC 
(see right panel in Fig. [2J. This may explain why the so urce is 
not specifically listed in Table B.l of lUeda et al.l (l2012h . There 
is no other CO peak in the SGMCs with an obvious counterpart 
in the near-IR, except for the embedded SSC in SGMC 4/5. 

Fig. [3] shows a comparison of the CO and CRIRES high- 
resolution (6 km s 1 ) H 2 1-0 S(l) spectra of PCC (Hll). The 
CO spectrum is the peak emission at the position of the PCC cor- 
rected for the surrounding emission, measured over an annulus 
outside the source, to isolate the CO velocity component asso- 
ciated with the PCC. CO(3-2) and H 2 1-0 S(l) spectra are re- 
markably similar. Table 1 lists the parameters of our Gaussian fit. 
The H2-luminous PCC has a velocity dispersion (40 km s -1 ) sig- 
nificantly higher than those of GMCs with the same sizes in th e 
Milky Way (5 km s^ lFalgarone et alJl2009l:lHever et al.ll2009h . 

Using the velocity dispersion of CO(3-2), we obtain a virial 
mass of M vir = 5Rcr 2 /G = 4.6 x 10 7 x (o-,,/40 kms -1 ) 2 M Q . 
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Table 1. CO(3-2) properties of the SGMCs (top) and PCC (bottom). We include the integrated H 2 1-0 S(l) fluxes, and the 
H2 1-0/2-1 S(l) and H2/CO flux ratios. For PCC, the line parameters are measured with a single-component Gaussian fit. 



Source 


velocity 
component 


Vlsr 
km s" 1 


km s" 1 




J CO 

Jy km s" 1 


M , 

M Q 


Fn 2 1-0 S(l) 

erg s" 1 cm -2 


Fu 2 2-1 S(l) 
1-0 S(l) 


H 2 1-0 S(l) 
CO(3— 2) 


SGMC 1 


a 


1417±2 


105+6 


380+24 


5.6+0.3 x 10 8 


1.6xl0" 14 


0.2 


2.1 


b 


1521+1 


63+3 


289+16 


4.2+0.2 x 10 s 


SGMC2 


a 


1469+1 


90+2 


289+8 


4.2+0.1 x 10 s 


8.3xl0" 15 


0.2 


1.6 


b 


1613+2 


114+5 


176+10 


2.6+0.1 x 10 8 


SGMC 4/5 


a 


1505+1 


58+2 


171+9 


2.5+0.1 x 10 8 


1.3xl0" 14 


0.3 


2.6 


b 


1586+3 


126+6 


256+12 


3.7+0.2 x 10 8 


SGMC 6 


a 


1506+2 


50+8 


48+12 


0.7+0.2 x 10 8 


6.7xl0" 16 


<0.5" 


0.2 


b 


1561+5 


113+6 


216+15 


3.2+0.2 x 10 8 


H 2 source (PCC) 




1534+2 


93+4 


7.6 ±0.4 


1.0+0.1 x 10 7 


7.6xl0" 16 


0.1 


8.7 



8 H2 2-1 S(l) flux corresponds to an upper limit estimated from the noise. 
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Fig. 3. (left) CO(3-2) emission from SGMC 2, where emission 
between 1350-1540 km s" 1 and between 1540-1750 km s" 1 
are shown in blue and red, respectively. Contours show the 
H 2 1-0 S(l) morphology, (right) CO(3-2) (red) and CRIRES 
H2 (blue) spectrum of PCC. The H2 spectrum is smoothed to 18 
km s" 1 resolution. The bar corresponds to 77 mJy beam 1 for 
CO(3-2) and 4.4 xlO" 17 erg s" 1 cirT 1 A" 1 forH 2 1-0 S(l). 



Since ALMA does not spatially resolve the PCC, we instead 
used the 50 pc size measured with SINFONI. The exceptionally 
high H2 1-0 S(l)-to-Bry line ratio (>15, Hll) provides unam- 
biguous ev idence that the H2 emission of the PCC is powered 
by shocks dPuxlev et al.l [l990). The H2/CO ratio, i.e. the energy 
dissipation rate per unit mass, is also exceptionally high, a factor 
5 higher than that of the SGMC 2 complex overall. 

PCC appears to be located at the interface between blue and 
redshifted gas (Fig. [3) where CO shows a steep velocity gradient 
(~1 km s" 1 pc" 1 in the position-velocity diagram in Fig. [2]). The 
observed properties of PCC are consistent with a scenario where 
the formation of SSCs is triggered by interactions between two 
gas flows. In SGMC 2, depending on the full three-dimensional 
geometry, the flows could either be colliding or creating a large 
velocity shear, and most likely a combination of both. In either 
case, the interaction drives a turbulent energy cascade in which 
kinetic energy is being dissipated. This is where we would ex- 
pect the highest energy dissipation rate. 

The bolometric luminosity of the PCC is ~10 7 L Q . 
Observations of, e.g., NGC 1333 (M aret et al.ll2009h show that 
the bolometric luminosity of protostellar outflows is on the or- 
der of 10 5 L Q x M w j n d, where M w j n d is the stellar mass loss rate 



in M Q yr" 1 . The small embedded stellar mass of M*=4 x 10 4 M Q 
(HI 1) makes protostellar winds an implausible energy source. 

The cloud luminosity may be accounted for by the dissipa- 
tion of the cloud kinetic energy for a cloud mass of a few 10 7 M Q 
- a value comparable to the virial mass 5 x 10 7 M - and a dissi- 
pation timescale of 1 Myr. This is comparable to the cloud cross- 
ing time, and also the dynamical time scale associated with the 
velocity gradient of 1 kms" 1 pc" 1 at the position of the cloud. 
The similarity of both timescales indicates that the cloud may 
still be forming by accreting gas, and therefore that a significant 
part of the cloud luminosity may be powered by gas accretion. In 
any case, the time during which the PCC is a bright H2 emitter 
is short, about 1 Myr. 

This short timescale may explain why we do not find more 
bright compact sources in our new H2 data. IWhitmore et alj 
( 20 1 0[) list five massive (>10 5 M ), young (< 5 Myr) SSCs over 
the part of the overlap region covered with both ALMA and 
SINFONI (Fig. [TJ. This is consistent with finding only a single 
bright PCC with the SINFONI data if the H 2 -luminous phase 
does not last longer than a few Myr. 

Our analysis gives a foretaste of the power of combining 
mass and energy tracers to study the dynamical state of molecu- 
lar gas in galaxy mergers and the early stages of the formation of 
SSCs. In the future, this approach can be extended with ALMA 
and VLT using additional tracers of energy dissipation and mass. 
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